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I. INTRODUCTION 


A. BACKGROUND 

The Captured Air Bubble surface effect ship presents 
a new concept in shipping for the surface force of the 
U.S. Navy. Within a few years these vessels, supported 
partially by a cushion of air and partially by the more 
conventional hull buoyancy forces, will be joining the 
Fleet. In preparation for that event it is of extreme 
importance for detailed studies to be undertaken so that 
the capabilities and limitations of such vessels may be 
Tul)yY understood: 

The term surface effect ship applies to a broad spec- 
trum of dirtferent types of vessels. Among them are the 
Air Cushion Vehicle (ACV), the Ground Effects Machine 
(GEM), the Semi-Submerged Ship (SSS), the Captured Air 
Bubble (CAB), and others. This study is concerned with 
the CAB surface effect ship. 

The Captured Air Bubble craft consists of two rigid 
Sidewalls which resemble catamaran hulls and flexible 
Seals between the water surface and hull structure across 
the bow and stern of the crata., The enclosed volume is a 
plenum chamber into which air is blown by a lift fan system. 
The pressure in the plenum chamber acting over the area 
of the plenum provides a lift force sufficient to support 


amlarze Proportion of the total crazt wergh>. The remainder 
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of the supporting force is made up by conventional hydro- 
dynamic forces acting on the sidewalls, plus some lift 
forces due to the bow and stern seals and aerodynamic 
lift at higher operating speeds. 

In the "on bubble" mode of operation, at the higher 
operating speeds, approximately two-thirds of the total 
craft weight is supported by the air cushion, and therefore 
relatively shallow sidewall immersion is required to support 
the remaining weight. The shallow draft produces a lower 
hydrodynamic drag and thus higher speed for a given power 
than could be obtained for a displacement hull. 

The XR~3 is a CAB surface effect ship of approximately 
three tons displacement. The XR-3 was built by the Naval 
Ship Research and Development Center (NSRDC) in 1966. 
Following periods of testing by NSRDC, Annapolis Division 
NSRDC, and Aerojet-General Corporation, the XR-3 was de- 
livered to the Naval Postgraduate School (NPS), Monterey, 
California in Marchal970 forthe anvestigation of basic 
and advanced surface effect ship technology. 

itiivadavoulon bo -unee.R=5 “tnerencxrse -Otner GAR test 
craft. Among these are the Aerojet-General 100-A and the 
Bell Aerospace Systems 100-B, both of approximately i100 
tons displacement. These two vessels were constructed as 
prototypes of a 2,000 ton ship. To date speeds in the 80- 
knot range have been attained by these vessels. 

Undervwecl.s Government contract, “Oceanircs Incorp— 


orated developed a digital computer simulation of the 
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CAB ship loads and motions for the Surface Effect Ship 
Project Office (SESPO) in Washington D.C. This loads and 
motions program with the mathematical model of the Bell 
Aerospace Systems 100-B with six-degrees-of-freedom has 
been in use at the W.R. Church Computer Center at NPS 
Since October 1972. 

In December 1973, LCDR. D.G. Leo, USN and LT. R. Boncal, 
USN produced a mathematical model with six-degrees-of- 
freedom of the XR-3 test craft derived from the above- 
mentioned simulation model of the 100-B (Ref. 1). The 
modular construction of the 1i00-B simulation program 
facilitated the re-modeling of the program where necessary 


to fit the design differences between the 100-B and the XR-3. 


B. OBJECTIVES 

HE 15% the ob jective of this thesis to: validate the 
non-linear, six-degrees-of-freedom, mathematical model of 
the XR-3 CAB test craft in a comprehensive manner by direct 
comparison between instrument-measured variables on the 
XR-3 and computed variables produced by the simulation 
peoframs. “in the process of carryins out-this objective, 
subroutine modifications were made to the simulation 
program to improve on the structure of the simulation model 
and justification is given in the sections that follow for 
the changes that were made. 

The data taken on the XR-3 measured in real-time, 


analog fashion using FM modulation was filtered, converted to 
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digital form, and processed in such a manner that it could 
be interfaced directly with the simulation program and a 
time history of true differences for each variable could 
be produced at each operating condition. 

The present study is limited to operation of the test 
craft in calm water conditions at the test site located 
at Lake San Antonio Reservoir, California. The range of 
operation of the test craft includes: 

1. Straight runs at various speeds with various 

locations of the center of gravity. The center 
of gravity was shifted by placement of two 125 
pound water tanks at different locations on the 
deck™oi the test crai tl. 

2. Straight runs at various speeds in an unladen 

COnd1 t10n. 
3. Runs made at various speeds with various rudder 


angles applied to validate turning performance. 
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Ii. COLLECTION OF DATA 


A. INSTRUMENTATION ON THE XR-3 

In order to support the continuing studies being under- 
taken at the Naval Postgraduate School on Captured Air 
Bubble vessels, the XR-3 has been instrumented with devices 
capable of measuring and recording time histories of many 
different variables for later, post-run analysis (Ref.2). 
The general configuration of the XR-3 is shown in Figure l. 
Table I shows a listing of the measurable parameters, the 
range of measurement, and the accuracy with which the 
measurements may be made. The locations of the measuring 
instruments are shown in Table II. The measured data are 
recorded on a Pemco Model 120-B magnetic tape recorder 
which features 14 tracks for recording data and an edge 
track for voice annotation during measurement runs. Each 
of the 14 tracks on the tape recorder is capable of record- 


ing = 1 volt RMS with a tolerance of i O57. 


B, @0ORER SUPPLIES 

The tape recorder is powered from a 110 volt 60 cycle 
auxiliary power unit, located on the aft weather deck as 
Siowa Ine Picure Pays throug sa Frencota.c.: O..d:C. Converter 
which provides a 26 volt d.c. power source. The three 
accelerometers for measuring vertical, lateral, and fore- 
and-aft accelerations are powered from a 12 volt storage 


Dattery through a d.c. to d.c. converter which provides a 
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Table I. 


Parameter 


Pore Jnrus v 
starboard Thrust 
Velocity 

Height 

Bowseal Press. 
Plenum Press. 
Sternseal Press. 
Pitch Angle 
Roll Angle 

Yaw Angle 

Pitch Rate 

Roll Rate 

Yaw Rate 

Rudder Position 


Longitudinal 
Acceleration 


Lateral 
Acceleration 


Vertical 
Acceleration 


*Height is used to compute Draft as 


XR-3 Measurable Parameters 


Simulation Program 


Name 
THS TP 
THS TS 
VEL 
None* 
PBARB 
PBAR 
PBARS 
THETAR 
DERL 
DPSI 
PDEG 
QDEG 
RDEG 


DELRS 


None 


ACCLAT 


ACCEL(3) 


18 


Range 


0-500 lbs 
0-500 lbs 
O-40 kts 
= 2 feet 
0-60 psf 
0-60 psf 
0-60 psf 
z 15 degs 
a5 15 degs 


= 480 degs 


730 degs/sec 
+30 degs/sec 


“30 degs/sec 


i+ 


45 degs 


0.2g units 


Oi22 Units 


Accuracy 


i+ 


eee 


t+ 


eve 


I-+ 


1 kt 


i+ 


Oo 


I+ 


> pst 


I+ 


0.5 pst 


‘ce 


O35. psi 


t+ 


0.5 deg 


I+ 


0.5 deg 


i+ 


0.5 deg 


“0.5 deg/sec 
BOs deg/sec 


to.5 deg/sec 


not available 
not available 


not available 


0.8 to +1.2g not available 


explained in text. 





Table II. 


sensor 


Pressure 
Transducers 


Bow Seal 


Plenum 
Stern Seal 


Velocity 


Height 


Gyros 
Pi cen 
Row 
Yaw 


Port Tarus 6 


Stardoara Thrust 


Vertical 
Accelerometer 


Lateral 
Accelerometer 


Longitudinal 
Accelerometer 


Distance Forward 
of Transom 


18" 
12'4s" 
1'6" 
Doe 


tod" down from 
deck 


2a 
ou down from 
deck 


1286 
is). 
Wo 
0 


0 


13 t 63" 


pox 6" 


13" 53" 


1) 


Location of Sensors 


Distance from 
Center Line 


2'hen 
athe" 
274s" 
i 6" 
Gitons 
1 *o2" 
1 oO" 


a0 t Bin 


POEG 
Port 
Port 
POR 


> tarboard 


Port 


ror. 


ror. 





28 volt d.c. supply. The three gyros which measure pitch, 


roll, yaw, pitch rate, roll rate and yaw rate are powered 


by the 12 volt storage battery. The gyros also receive a 


1 volt d.c. excitation signal from the 110 volt a.c. power 


supply fed through an a.c. to d.c. converter. 


The height sensor mounted forward of the bow is 


powered directly from the 110 volt a.c. auxiliary power 


unit. The height measurement signal is converted to a 


71.0 volt d.c. Signal by a height sensor conditioner. The 


velocity probe, 
sensor, extends 
testcraft. The 
PROvCen tS. Une wl2 


converter. The 


mounted on the same assembly as the height 
into the undisturbed water forward of the 
28 volt d.c. power source for the velocity 
vol. Datbwery*tnroush: ~ne d.c. oOo d.c, 


rudder position measurement is made by means 


Ot vas OaG Ce ll twitch receives “acl2-volLt dC. supply irom 


the battery. 


The two load cells for measuring port and starboard 


thrust and the three pressure transducers for measuring 


bow seal, plenum, and stern seal pressures are all powered 


Dy ae convo oso ee sUppd yarrom the 2 volt duc. Storage 


battery through 


une UGC. vO ate). S6onverter;, All ‘five 


UMitena | SsOerece1 6) a05 VO UV JC eC tatuon sienad Tron 


an amplifier-transducer package. 
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ITI. DATA PROCESSING 


A. ANALOG TO DIGITAL CONVERSION 

To provide measurements for validation of the mathe- 
matical model, data-gathering runs were made with the XR-3 
testcraft on four different days at the Lake San Antonio, 
California test site. The data was recorded on magnetic 
tape using FM modulation, then processed into a suitable 
form for comparison with corresponding values computed in 
the XR-3 simulation program on the IBM 360/67 Digital 
Computer. 

The data was converted from analog to digital form 
on an XDS-9300/CI-5000 hybrid computer at the NPS 
Computer Laboratory. From previous experience in data 
reduction by the XR-3 Research Group in the Aeronautical 
Engineering Department at NPS, it was found that filtering 
of the analog signal was required to remove unwanted high 
frequency noise from the recorded data. To provide this 
required filtering six Krohn-Hite digitally-tuned filters 
were assembled so that six tracks of recorded data could 


be processed simultaneously. 


B. DESCRIPTION OF EQUI BMENT 

The analog tape containing raw data from XR-3 test 
runs was mounted on a Honeywell Tape Recorder/Reproducer, 
Model 96 for playback. Connections were made through 


trunk lines to lead the output from the six selected tracks 
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to the COMCOR CI-5000 Analog Computer as shown in Figure 2. 
Each channel of recorded data was then passed through an 
amplifier with a gain of 5 before being led out of the 
CI-5000 to one of the filters. The filters were operated 
in the low pass-maximally flat mode which provided a four- 
pole Butterworth response with an attenuation of 24db per 
Sctave irom the cut-off frequency. Cut off frequencies 
were set at 10 Hertz for all channels except port and 
starboard thrust which were set at 1.0 Hertz. Each 

Signal at the output of its respective filter was led back 
into the Ci-5000 for another amplification by a gain factor 
of ten. A Brush 8-Channel Recorder was also patched into 
the system at this point so that each of the six selected 
data tracks could be continuously monitored. Figure 2 
shows the signal flow patch through the hybrid computer 
system for analog-to-digital conversion. 

For analog-to-digital conversion, digital-to-analog 
conversion, and other applications, the CI-5000 Analog 
Computer and the XDS-9300 Digital Computer are intercon- 
nected through trunks. For A-to-D conversion in this 
application a sampling frequency of 200 Hz was chosen. In 
steady-state calm-water runs with the XR-3 none of the 
measured parameters can be expected to vary at frequencies 
Sreater “wien 10 Hz. Thereiore; ‘selection of a 200 Hz. 
sampling frequency more than satisfies the requirements 


of the Sampling Theorem. 


£2 





NOTSUSANOD WIESTC GL CNW ONTand Old wvNIIS  °¢ JunoIs 


MOTS TWNOIS NOISYSANOD V/O SLVOIONI S3NI1 Q3HSYO 
MOTS TWNOTS NOISYSANO9 O/W 3LVOIONI S3NI7 G110S .:310N 






ddVl OVW 








> & YAddOIS4 
HSNYS 
VANNVHD 8 

















YALYIANOD 2 
¥ Ol Qa Serene > 
JLIH-NHOUD 
YALYIANOD 
aol Y 
YILNGWOID IWLIDIG Y4LNdNOD DOTWNY 
00f6 sax | 0006-19 YOINOD 96 GOH! YIINGONdaY 


TIS MASH OH 





The first step in the A-to-D process was to zero all 
possible bias errors in the system. Calibration steps 
are recorded on separate channel strips on the analog 
data tape which indicate zero and maximum signal readings 
for each sensor. Bias potentiometers at the input of the 
first stage amplifier on the CI-5000 were adjusted to 
properly zero all channels during the pre-test run cali- 
bration of the XR-3 sensors. Each signal was monitored 
on a digital voltmeter during the calibration procedure. 

Execution of the FORTRAN program listed in Appendix A 
on the hybrid computer system enabled selective conversion 
iFom analoe to-diegital form of-specifie time periods of 
analog data. Continuous monitoring of the voice edge track 
on the data tape vrovided a reference as to what the XR-3 
conditions were at all times. After each file of data had 
been digitized and written on the digital magnetic tape, 
the tape was re-wound and that file most recently recorded 
was run through the digital-tc-analog conversion. The 
D-to-A signal was monitored on the Brush Recorder and com- 
pared with the analog signal to verify the validity of the 
COnVersion process... “[neetinisned product-of tne. A=to—-D 
conversion was a reel of 7-track digital tape consisting 
of many files of digital data. Each file on the tape 
contained a time history of six selected variables measured 


Lon i speeiiic Opera tine hCONGiLLoOnson “tne! 2R= 3. 
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C. INTERFACE WITH IBM SYSTEM 360/67 

Introduction of the 7-track tapes into IBM System 
360/67 at the W.R. Church Computer Center at NPS entailed 
some further processing of data. In the XDS 9300 a binary 
number is represented by considering its octal representa- 
tion, made up of 3-bit digits, whereas in the IBM 360 
the same number is represented by hexadecimal, made up 
of 4-bit digits. As an example of this difference, a 


decimal number (668), 7S (1234) octal and (29C)4¢ 


( 
hexadecimal. It was necessary, then, to write a FORTRAN 
program to read each record ina file from the 7-track 
tape, make the conversion from octal to hexadecimal, store 
the data in hexadecimal form on a temporary storage disk, 
ange tom priit Out a DOortloneo: iach record in decimal form: 
A listing of the program is included as Appendix B. The 
program uses the pre-compiled assembler language subroutine 
FORM to convert the data to its nexadecimal representation. 
Figure 3 shows the flow of information through the IBM 
System 360/67 computer in this validation. 

The printed output from the conversion program is re- 
presentative of the data in its original measurement form 
multiplied by a gain factor of 50 in the A-to-D conversion 
process. At this point in the process the numbers printed 
out were compared with the signals recorded on the Brush 
Recorder during A-to-D conversion to insure that the 
accuracy of the data had been preserved. This was indeed 


the case with the exception of some small amplitude, 
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high-frequency fluctuations in the digitally reproduced 
data which did not appear to be present on the Brush 
Recorder tapes. 

The final phase in data processing was to recover the 
measured data into the units which correspond with the 
measurement, i.e. pounds of thrust, degrees of pitch angle, 
knots of velocity, etc. to facilitate direct comparison 
in the XR-3 simulation model. To accomplish the recovery 
of data it was necessary to read the converted data from 
the temporary storage disk, multiply each set of variables 
by the inverse of the appropriate scale factor used when 
recording the data on the testcraft, divide out the gain 
factor wor 5O/introduced an the A= to-D process, and store 
the result on a temporary disk for later access by the 
XR-3 simulation program. 

The height measurement made on the XR-3 during data 
collection runs does not have an immediately available 
COuUnLer part, a1). The: SR=Jeproeram as. iidicated an Table 1. 
ho vyery Close approximation. tovdratt can be: obtained from 
the height measurement through suitable scaling, however. 
Thea Galeulavvon-ofvdralt 17 the xkR=34s mulation procram 
ple smade wath Ppeterence to dratt atthe center of gravity 
of the craft. The height measurement is made at a sensor 
lG@cated 191” forward of, the reference center of gravity 
and 6.375" dom from the deck level of the test craft 
(see Figure 1). The measurement of the height of the 


sensor above the surface of the water will vary as the 


2/ 





craft pitches up and down at the bow. Since pitch angles 
are, for the normal operating range, less than five degrees, 
the small angle approximation may be made for the sine of 
pitch angle. The draft measurement 1s synthesized by 
solution of the following relationship: 

DRAFT = SH + Z - HT + D*SING 


Water height above keel at the center of 
gravity 


where DRAFT 


SH = Height of sensor above CG = 7.145" 
Ze= Distance. 61. CG above Neel = 30.46" 
HT = Measured height of sensor above water surface 
D = Distance of height sensor forward of CG = 191" 
D*SIN@ = Variation of measured height due to pitch angle, 
A FORTRAN program was written to accomplish the process 
of re-scaling to measured units and to provide the ciose 
approximation to a draft measurement. A time index was also 
included in the program to facilitate plotting each variable 
versus time on a CALCOMP plotter. This feature proved to be 
valuable as a final step in the process because each variable 
could be closely examined in graphical form prior to acces- 
Sing “tne data by the: AR=>) Simulation program: 12t was 2¢ 
this point in the procedure that some difficulties were 
encountered in early attempts at signal processing, as 


described in the followings section. 


D. FREQUENCY COMPONENTS ENCOUNTERED IN DATA 
Examination of the CALCOMF plots of the digitally 


reproduced measurement data revealed significant small 
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amplitude, high frequency noise impressed on the measure2- 
ment signal. The noisy signal was not immediately evident 
on the Brush Recorder tape, however on close inspection it 
was apparent that some noise was passing through the filters. 
Figures 4,5,6, and 7 are representative CALCOMP plots of 
processed measurement data of thrust, pitch angle, velocity 
and plenum pressure respectively. These graphs are in- 
cluded to show the effect of noise on the processed data. 
The source of noise present in the measurement data 
has not been established and can only be rationalized to 
be a combination of several possible environmental factors. 
First of these is vibration aboard the XR-3 testcraft from 
the twin engine propulsion system and fan drive system. 
Second, the testcraft 1s rarely run in completely glass- 
calm water conditions at the test site. Surface disturbances 
are nearly always present due to numerous pleasure craft, 
Wind-induced ripples on the lake, and occasionally, the 
Wake Of. the tes terattcatselt. “Suchdisturbances can have 
a pronounced effect on the motion of the XR-3. Finally, 
60-cycle power supply interference can be assumed to be 
present due to the antenna effect of unshielded wiring 


DOUNvdn Lie bes tCral tcand Im the NPS. Computer Laboratory. 


i. Elimination of Unwanted High Frequencies 
Recognizing that the noisy data rendered it all but 


unusable for any meaningful attempt at validation, methods 


were sought to eliminate the noise and to obtain a "clean" 
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Figure 5. Plot of Pitch Angle Vs. Time 
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Signal and, at the same time, to preserve the validity of 
the actual measurement. The first effort was to filter 
at a lower frequency. The result was displeasing for, 
while the magnitude of the impressed noise was somewhat 
reduced, lower frequency oscillations which represented 
actual craft motions were nearly eliminated. The next 
possible solution was to employ two filters in cascade, 
thereby obtaining a low-pass filter with 48 db per octave 
attenuation above the cut off frequency. This procedure 
accomplished the main objective of removing virtually all 
the noise and passing only signals representative of the 
craft motions. The disadvantage of this method was that 
only three instead of six measurements were now possible 
because of the hardware limitations and synchronization 
became an insurmountable problem. 

Me third possibility of “removine the noise from 
the measurement signal was to proceed on the assumption 
that the noise was gaussian and white and to employ a 
smoothing or time-averaging technique. The technique em- 
ployed was to process the data as originally set forth and 
then,;-1n the final scaling process on. the: digital computer, 
to obtain a moving time average. The algorithm implemented 
on the digital computer was to calculate the sum of the 
first ten points in the time history of a measured para- 
meter, divide the sum by ten and call the resultant figure 
the first point in the time history. Mext, a counter was 


incremented by five and the process was repeated to obtain 
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the next point. In this manner the data is reduced from 

200 samples per second to 20 samples per second. MThe over- 
lapping smoothing technique tends to preserve low frequencies. 
By comparison there were insignificant differences (less 

than two per cent) between this smoothing process and the 
aforementioned cascade filtering. The advantage realized 

was a more economical use of available equipment (filters) 

and a great savings in time required for data processing. 

A listing of the program used for re-scaling the 
data to its original measurement form and including the 
smoothing routine is included as Appendix C. Figures 8, 

9, 10 and il are CALCOMP plots of smoothed measurement data 
of thrust, pitch angle, velocity and plenum pressure 
respectively. By comparison with the measurement data 
plotted in Figures 4,5,6, and 7, it is evident that the 
moving time average smoothing techniaue has preserved the 
measurement information of value and has eliminated the 


impressed noise. 
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IV. INTERFACE WITH XR-3 SIMULATION PROGRAM 


A. LOGIC ADDITIONS TO SIMULATION PROGRAM 

The XR-3 Loads and Motions simulation program was 
modified to the extent that, on execution of the program, 
the processed measurement data could be accessed if desired 
and on-line computation of differences could be made. ‘The 
output from the program would then be a time history of 
the differences between computed and measured parameters. 
This entailed the addition of another option switch in the 
input data which is rea by subroutine INCON. If the option 
switch is on or equal to one, Subroutine INCON reads the 
measurement data into core from the disk where the data 
ISuStOred. Tf “thenOpulOneeWwltcn, 1s-Off Or ECqualk to Zero, 
no reading of measurement data occurs. 

It was decided to restrict validation or data compari- 
son runs to ten seconds of simulation time because of core 
storage and computer time allocations. This decision set 
the size of the measurement data arrays at 200 elements. 

A COMMON block was added in subroutine INCON and RHS so 
that the measurement data could be available in subroutine 
RHS. As with other COMMON blocks, the feature to zero all 
elements of the added COMMON block was included in the 
main program. The following ZAR-3 test craft variables 

and the associated simulation names were selected for 


Va bidation: 
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Totabt Thrus G ES 


Velocity VEL 
Height DRFT 
Pitch Angle THE TAR 
Roll Angle DPHI 
Plenum Pressure PBAR 
Yaw Angle DEPST 
Rudder Angle DELRS 
Vertical Acceleration ACCEL(3) 
Yaw Rate RDEG 


The necessary logic was added in subroutine RHS such 
that the appropriate element with respect to time in the 
measurement array could be compared with the computed 
variable and a difference formed. The resultant differ- 
ences were written into a temporary file for output by 
subroutine COLFIL after the simulation run reached the 
specified final time of ten seconds. Computation of 
differences was arranged such that positive values indi- 
cated that the computed parameter was greater than the 
measured parameter for any given sample time. The print 
interval for output purposes was set at one point every 


0.05 seconds to correspond with the measurement data rate. 


Be Me tHOD OF THRUST INPUT AND OSSERVED: DIFFERENCES 
The first step in the validation process was to 
initialize the mathematical model at a particular operating 


point with the expected trim conditions of velocity, pitch 


4 





angle, draft, and plenum pressure. The feature exists in 
the simulation program to provide a thrust map of up to 

25 points. Additional logic was included to select every 
fourth point in the 200 element array of measurement 
thrust to be applied to the mathematical model as an input 
variable and a ten second simulation proceeded to produce 
computed velocity. All other measurement parameters except 
velocity were ignored during initialization runs. The 
difference between computed velocity and measured velocity 
was taken as an output from the model. Following initial- 
ization of the model, another computer simulation was made 
using computed thrust as necessary to hold velocity con- 
stant. All measurement data for the run was read in and 
differences between computed and measured parameters were 
taken as output from the computer model. 

For validation studies using measurement data with 
Very Nearly constant Thrust applied? on the testcrait, the 
henesoi1ns procedunmeycouldys be-shnor tened "by zmcluding only 
the last phase. Provided that the initial trim conditions 
were knovm for input to the computer model, no loss of 


accuracy resulted an the. process of validation. 


Cee THOL eCr SRUDDER WAND 2takUs il INPUD. VO TE MOREL 
Similar to the thrust map provision is a provision 

fOr a) rudder map anput to. the computer model... “Up to 25 

points may be included in a rudder map and the necessary 


logic was included in the program to select every fourth 
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point in the 200 element array of measured rudder angle 

for input to the mathematical model. The objective here 
was to apply the same inputs to the mathematical model as 
were used on the testcraft along with the initial trim 
conditions for each case. This method proved to be of 
little real value because a steering cable failure to the 
port propulsion engine on the testcraft rendered the rudder 
position measurement unreliable. To substitute for direct 
steering control on the port engine, the starboard engine 
was mechanically linked to the port engine by means of a 
long, rigid bar. Thus, the starboard engine controls the 
port engine in the rudder position and alignment of thrust 
VECtors: On thestestcrait. Virtually ali precision in rudder 
angle measurement is lost due to unavoidable play in the 
linkage and lack of direct Suserime control Over both 
engines. 

Validation studies involving turning maneuvers with 
the 4R-3 testcraft and the matnematical model were carried 
out by matching yaw angles and yaw rates for a spectrum of 
operating conditions and by noting the differences between 
computed and measured rudder angles for each operating 


Cond ftom; 
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V. PROGRAM SUBROUTINE MODIFICATIONS 


Aa REPLACEMENT OF “CENTER VOF- PRESSURE CURVE 

Newman and Poole reported in Ref. 3 the existence of 
a pressure wave and the attendant gradual deformation of 
the surface of the water due to a moving pressure distribu- 
t10n With regard to CAB type hulls, NSRDC found that, in 
towing tank tests of CAB scale models, the craft tends to 
pivot in pitch motion about a point approximately 75-80% 
of the wetted length of the sidewall forward of the tran- 
som (Ref. 4). Based on these findings, it has been 
postulated that the center of volume within the plenum 
chamber moves aft with increasing speed. 

Accordingly, the calculations made in the XR-3 sim- 
Wlatlon profran-to account for a shitting center of pressure 
were removed from subroutines RHS and WAVES. Computations 
were inserted in the program to account for the addition 
of an incremental wedge of volume to plenum volume which 
varies with craft speed. What follows is a description 
by subroutine of changes which were entered and the basis 
for the inclusion of the changes. 

1. Subroutine INCON 

A previously unused space on the input data card 
number 00701 was used to read in the location of the pres- 
sure Wave pivot point, measured in feet and tenths of 


feet forward of the transom. The statement XPWV = TEMP(3) 


tpl 





was inserted in Block 07 in the INCON subroutine. A value 
of 17.2 ft. was used for XPWV in this study. 

The relationship for wave drag produced by the move- 
ment of the pressure bubble over the water surface is, 
from Ref. 5:. 


Anubble = - (P, - P) ¢ 


The same formula is mechanized in the simulation program as: 
FXPWAV = -PBBAR * CF * PWVCON 
where CF = 0.37 / ((U/FNCON) **1.5655981) and PWVCON = 
4.0% WEIGHT /(RHO *G* XLBW). The coefficient 0.37 and the 
exponent 1.5655981 are correct for the length to beam ratio 
of the XR-3 as determined from Ref. 3. The variable name 
WATSLP is given to the slope of the water surface inside the 
plenum and is calculated as: 
WATSLP = PBBAR * CF * PWVCON / WEIGHT. 
WATSLP represents the angle at the apex of the wedge-shaped 
increment of volume to be added to plenum volume. Since 
the ratio of PBBAR*CF*PWVCON ( i.e. bubble drag) to WEIGHT 
will always result in a small number for speeds above hump 
speed, the small angle approximation may be used without 
loss of accuracy. Therefore, the volume increment due to 
the water slope to be added to plenum volume is given by 
C2 oe eA OLE i Kea Ae 
where XL is the length of the plenum chamber and AB is the 
area of the plenum at the water surface. The quantity 0.5 
« WAISLP * XL is then the vertical cross sectional area 


of the wedge-shaped volume increment. 
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Two other calculations were added in subroutine 
INCON for use in other subroutines during the simulation. 
The distance from XPWV from the forward extremity of the 
plenum is calculated as 
XLXPWV = XLBW - XPWV 
where XLBW represents the plenum length at the water 
surface. Also the distance between the craft center of 
gravity and XPWV is calculated as 
XPWVXS = XPWV - XS 
where XS is the longitudinal center of gravity measured in 
feet forward of the transom. 
2. Subroutine BOWSL 
In the computation for ELSKI, the wetted lenght of 
the bow seal, an incremental length is added to the wetted 
length of the seal. The increase is due to the water slope 
and 1s calculated simply as 
XLAPWV * WATSLP 
and added to the existing expression for ELSKI. 
o> * SubrOUTINeE SSTNSL 
Similar to the wetted length correction in BOWSL, 
an incremental correction is computed and applied in 
Subroutine sINol. Here the correction. 1c opposite in 
Slen and is applied to the wetted length of the stern seal, 
ELSKI as 
-XPWY * WATSLP. 
4, Subroutine RHS 
The water slope correction is computed in subroutine 


RHS as a continuous variable over the period of any given 
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Simulation time. Recognizing that WATSLP is calculated as 
the ratio of bubble wave drag to total craft weight and 
that bubble wave drag varies with velocity and plenum 
pressure, WATSLP must be considered as a continuous 
variable. Accordingly, the following computation was in- 
teoduced in Ris: 

WATSLP = -FXPWAV / WEIGHT 
The plenum volume is corrected for the water slope on each 
pass through subroutine RHS by the addition of a volume 
increment as in subroutine INCON: 

Peso WAS LE a he AB . 

5. Subroutine SIDEWL 
The addition of a wedge-shaped volume to plenum 

volume required an alteration to the computation of the 
drafts along the sidewalls inside the plenum chamber. 
Displacement of water is calculated as 

PBHEAD = PBAR / (RHO * G) 
where PBAR is the instantaneous plenum pressure. The im- 
mersion depth on the inside of each sidewall section is 
Caleulated 17 ethe relation 

DDIN = DD - WATSLP * (XPWVXS - XX(,K)) 
where DD represents the uncorrected water level height 
above the keel along each sidewall section and XX(J,K) 
me thevdis tance trom tne’ transom of the center,of cach «of 
the eleven stations along the sidewall. 

A revision of the method of computing sidewall gaps 


at any given station was also necessary. A logical IF 
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statement was included to test each value of DDIN after 
it was computed. If DDIN is less than zero another logical 
IF statement must be tested to determine if a gap exists 
between the bottom of the sidewall (keel) and the water 
surface. If DSW (J,K), or the actual sidewall immersion 
depth, is less than PBHEAD, then a gap for that sidewall 
section is computed as 

GAP (J,K) = - DDIN * (1.-(DSW(J,K))/ PBHEAD) 
where the computed GAP for any given station is in units 
of feet. 

Further background information of use in this 


study is contained in References 6-11. 


Boe TER JOP VECTORS ALN: SUBROUTING -ROP 

Tne propulsion system in the XR-3 consists of two 
forty horsepower outboard motors. The port screw has a 
right hand rotation and the starboard screw has a left 
hand rotation. Steering control is obtained by rotation of 
the motors in the horizontal plane. Thus, there is some 
rudder effect from the foil-shaped lower motor housings, 
however the predominant turning moment results from re- 
alignment of the propulsion thrust vectors when rudder angle 
1s applied. In addition to the main thrust vector, directed 
ie tie wLiverot oral: cumouron. tor zero, rudder angle, there 
exists a side thrust vector for each propeller. Since the 
propellers are counter-rotating, and due/to- the direction 


Of FOtaAGiGN Of Cach propeller, the side’ thrust vectors are 
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oppositely directed and, for zero rudder angle, are 
directed inboard. 

The assumption is made that both port and starboard 
propulsion motors rotate together through the same rudder 
angle. It is also assumed that the magnitude of the main 
and side thrust vectors are independent of rudder angle. 

The propeller thrust diagram shown in Figure i2 was 
used to verify and correct the force and moment equations 
for the XR-3 propulsion system. The figure displays the 
Main thrust and side thrust vectors for both propellers 
Postilonedswitimrienty (posislVe)-ridder angles, positive 
pitch angleé, and negative (port side down) roll angled. 
The thrust vectors act through the port and starboard pro- 
Beller postions, Ars Py sand 2b trem tne crait center of 
gravity in the Body Reference coordinate system. 

The forces associated with the propeller alignment of 
Figure 1 aresshovwm vectorially in Figure13 in the Local 
Level coordinate system after having been corrected for 
Ditenm and roll@anscies. Whe assumpvul0on 1s made that pitcn 
angle and roll angle are limited to small angles such that 
siné=@ and cos@= 1, sind =@® and cos @ = 1, The follow- 
ing force equations result: 

he oS Eee Cos 
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PP2Y = THST2 * sin § 


STLY = -STHST1 * cos § 
ST2Y = STHST2 * cos § 
T1Z = THST1 * cos& * 9 
T2Z = THST2 * cos& #9 
STLZ = STHST1 * cos§ * @ 
ST2Z = STHST2 * cos$ *¢ 


Summing the forces in eacn of the three corrdinate 


directions yields the following result: 


FXS 
EAP 
1a gs 
in ay 
FZS 
Ray 


The 
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total force equations are found by summation of 


the force equations for the port ard starboard propulsion 


systems in each of the three coordinate directions as 


follows 
FX 
FY 
FZ 


EAP SEAS 
| Ob sees sae on 


BZP = Zs 


Corrections were made 
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Vi. DISCUSSION SAND EVALUATION OF RESULTS 


A. OBTAINING AGREEMENT AMONG VARIABLES 

Computer simulation runs in the initial stages of this 
study revealed that plenum pressure was consistently ap- 
proximately 2.5 psf higher than that recorded on the KR-3 
testcraft. Among other variables, thrust on the testcraft 
was higher than the computer model indicated for any given 
speed and calculated draft was greater by approximately 
three inches than the computer model produced. The decision 
Was made to first obtain agreement in plenum pressure 
between the mathematical model and the testcraft. 

The bow seal leakage area input to the mathematical 
model was listed in Ref. i as one of the parameters not 
known to be exact. Tne estimate of leakage area was re- 
ViScdrironO-0C.Ssaqsnt us tOrd.,10 sq. ft4 This slight alter= 
ation had the desired effect of reducing plenum pressure 
to agree with that measured on the XR-3 testcraft. An 
unexpected but welcome additional effect was nearly perfect 
agreement in thrust and draft. Aside from the alterations 
to the structure of the mathematical model there were no 
further changes made to the model during the course of 
the present study. 

i. Se toerent Rums Made With Ditierent €.G...Posi tions 

Validation runs were made interfacing data measured 


on the KR-=3 testcraft with the computer model... ‘The data 
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presented in Table III is a summary of average differences 
calculated between time histories of actual craft motions 
and simulated craft motions for the thrust, velocity, 

pitch angle, plenum pressure and draft variables. Also 
included is percentage differences for ease in qualitative 
comparison. It should be noted that percentage differences 
in pitch angle measurement appear to be quite large, how- 
ever in most cases the difference is within the measurement 
accuracy limitation of " 0.5 degree aboard the testcraft. 

The total craft weight on the measurement runs 
Micced ne Table. lll was 0270: lbs, Invaddition to the pilot 
there was a passenger in the co-pilot position and 250 lbs. 
of water in plastic containers aboard used for changing 
the location of the center of gravity of the XR-3. 

The data presented in Table IV is a summary of 
average time history differences with a fixed center of 
gravity on the XR-3. On the measurement runs listed in 
Table IV the total craft weight was 6050 lbs. Vertical 
acceleration was substituted for draft in this series of 
validation runs. 

Pieures is oe VO I 7.sand fo are typical of the 
graphical output from the XR-3 simulation program arranged 
to display tne time histories of differences between com- 
puted and measured variables for this study. 

2. ‘Rune Bade yo th Tunas 

A series of eleven validation runs were made with 


various rudder angles at various speeds. <A summary of the 
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average differences calculated between actual craft motions 
and simulated craft moticns during the turning maneuvers 
for thrust, velocity, pitch angle, roll angle, yaw angle, 
yaw rate and rudder position variables is presented in 
Table V. 

The relatively large differences in rudder angles 
in many runs is attributed to the lack of confidence in 
an accurate rudder position measurement. in the future 
Wiel drrect=s Tecrine contol, G1 both propulsion motors. is 
restored, it 1s recommended that confirmation be made of 
the results presented here with regard to rudder position 
and other variables in turning maneuvers. 

Bicures 19. 20; 2btand 22eare typreal of the <rapn- 
ical output from the XR-3 simulation program which display 
the time histories of differences in roll angle, yaw angle, 


yaw rate, and rudder angle, respectively. 


B. LIMITED TO CALM WATER CASES 

It 1s recognized that much of the measurement data 
Obtatnedeat tie An—>) Lest -siterat lare San Antonio, Cali- 
fornia contains small perturbations due to surface 
disturbances on the lake. No attempt has been made to 
measure the surface disturbances or their effect on the 
testcraft. Validation of the mathematical model in the 
presence of waves of any amplitude and wavelength is sug- 


gested as an area for further study. 
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Figure 19. Plot of Roll Angle Difference Vs. Time 
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A listing of the FORTRAN program for simulation of 
the XR-3 is included as Appendix D. A sample input data 


deck is also included. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


This study has shown that the nonlinear mathematical 
model with six degrees of freedom adequately describes the 
motions of the XR-3 in calm water conditions. The variables 
used in the validation included thrust, longitudinal velo- 
city, pitch angle, plenum pressure, draft and vertical 
acceleration in straight runs. For validation of turning 
dynamics the variables included thrust, longitudinal velo- 
city, pitch angle, roll angle, yaw angle, yaw rate and 
rudder angle. Differences that were found to exist between 
the computer model and the testcraft were less than 5% for 
the majority of the variables compared and those differ- 
ences that exceeded 5% are considered acceptable, given 
the range of the measured variables and the accuracy of 
the sensors aboard the testcraft. 

During the course of this study, it became apparent 
that the computer model exhibits a more oscillatory, or 
underdamped, response in pitch and roll motions than does 
ice Anaya ves LCA bo? AMCOncurrent study =by LCDR. Ruka Painiey, 
USN Upon] this Aspect icr “the= simulatiom with attention: to 
the modelling of the bow seal and the lift fan system 
25 Peported in «Reference 12. 

Hae method Of data processing setciorth. ine this study 
to produce data in adequate form for interface with the 


mathematical model is recommended for use in future 


(ea 





validation efforts. A considerable amount of time and 
effort was spent in signal analysis and data processing 

to prepare the measured data in a form useful for compari- 
son with the mathematical model variables. In order to 
eliminate the requirement for smoothing the re-scaled 
data, it is recommended that a study be made using a lower 
sampling rate for the A-to-D conversion to see if the high 
frequency noise can be eliminated from the digitally re- 
produced data. It is also recommended that a vibration 
study be conducted aboard the XR-3 in an attempt to deter- 
mine if the source or sources of the high frequency noise 
is the result of internal ship vibration acting on the 
measurement transducers, gyros, and accelerometers. fThe 
current data acquisition system installed in the testcraft 
should be adequate for such a study. 

A very fertile area for future study would be valida- 
tion of the computer model in the presence of waves. 
Instrumentation exists to measure wave conditions at the 
XR-3 test site during measurement runs and the simulated 
WaveSecCOULd DeMInpUt Toethne mathematical model. Direct 
comparison of craft behavior and simulation model calcula- 
tions could then be obtained to produce an evaluation of 


the accuracy of the simulation in sea states. 
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APPENDIX A 


LISTING OF FORTRAN PROGRAM FUR ANALOG TO 


DATA 


DIGITAL CONVERSION OF MEASURMENT 
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APPENDIX B 
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